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Abstract: Models of dark matter with ~ GeV scale force mediators provide attractive 
explanations of many high energy anomalies, including PAMELA, ATIC, and the WMAP 
haze. At the same time, by exploiting the ~ MeV scale excited states that are automatically 
present in such theories, these models naturally explain the DAM A/LIBRA and INTEGRAL 
signals through the inelastic dark matter (iDM) and exciting dark matter (XDM) scenarios, 
respectively. Interestingly, with only weak kinetic mixing to hypercharge to mediate decays, 
the lifetime of excited states with 5 < 2m e i ec tron is longer than the age of the universe. 
The fractional relic abundance of these excited states depends on the temperature of kinetic 
decoupling, but can be appreciable. There could easily be other mechanisms for rapid decay, 
but the consequences of such long-lived states are intriguing. We find that CDMS constrains 
the fractional relic population of ~ 100 keV states to be |o 10~ 2 , for a 1 TeV WIMP with 
a n = 10 -40 cm 2 . Upcoming searches at CDMS, as well as xenon, silicon, and argon targets, 
can push this limit significantly lower. We also consider the possibility that the DAMA 
excitation occurs from a metastable state into the XDM state, which decays via e + e~ emission, 
which allows lighter states to explain the INTEGRAL signal due to the small kinetic energies 
required. Such models yield dramatic signals from down-scattering, with spectra peaking 
at high energies, sometimes as high as ~ 1 MeV, well outside the usual search windows. 
Such signals would be visible at future Ar and Si experiments, and may be visible at Ge and 
Xe experiments, although 7-rays associated with nuclear excitations would complicate the 
signal for these heavier targets. We also consider other XDM models involving ~ 500 keV 
metastable states, and find they can allow lighter WIMPs to explain INTEGRAL as well. 



1. Introduction 

The range of ideas for dark matter has recently expanded dramatically. Rather than limiting 
ourselves to candidates such as the MSSM neutralino, there has been a profusion of models 
with new forces, new annihilation modes, and multiple states. The motivation for these 
new properties has come from a variety of sources 1 , but dominantly from dramatic new 
astrophysical signatures, among them, PAMELA [2], WMAP [3], ATIC [4] and PPB-BETS 
[5]. 

Perhaps most striking among the recent signals is the PAMELA positron excess [2] . The 
sharp rise of the positron fraction above 10 GeV seems inconsistent with interactions of CR 
protons with the interstellar medium, requiring a new primary source of positrons. This 
confirms the hints seen at HEAT [6, 7] and AMS-01 [8]. However, interpreting these signals 
as coming from dark matter annihilation is challenging in the traditional thermal WIMP 
scenarios. Annihilations to light leptons are often helicity suppressed [9], while annihilation 
to hadronic or gauge boson modes provides spectra that are generally too soft to fit the data 
[10, 11] 2 . Moreover, the copious anti-protons which are produced in these annihilations are 
generally an order of magnitude [11, 13] above the level seen by PAMELA [14]. Finally, the 
cross section required to explain the positron excess is typically an order of magnitude (or 
more) larger than the thermal cross section of 3 x 10~ 26 cm 3 s _1 . 

A simple solution to all of these problems arises from postulating a new light force for 
dark matter [1]. Annihilations of dark matter into light bosons which decay to electrons and 
muons can easily yield a hard lepton spectrum without significant anti-protons or 7r°'s [15] and 
provides a good fit to the PAMELA data [16, 17, 18, 19]. Moreover, the light mediator can 
provide an enhanced cross section at low (v ~ 220km/s) velocities through the Sommerfeld 
effect [1] or a capture into WIMPonium [20] 3 . 

Such a model can also explains other astrophysical anomalies. Specifically, such annihila- 
tions can explain [15, 25] the WMAP "haze" [26, 27, 28], understood as high energy positrons 
and electrons from WIMP annihilation, synchrotron radiating in the galactic magnetic field 
[27, 29]. Excesses of electrons and/or positrons in the 400-700 GeV range reported by ATIC 
[4] and PPB-BETS [5] are also naturally explained [25, 18, 19] by the same annihilation 
channels leading to the positron excess in PAMELA. 

If this new force is a gauge force, a very important consequence is the presence of quasi- 
degenerate states of dark matter [1]. While such states can have important consequences 
for direct detection signals [30] or cosmic ray production [31], it is important to recognize 
that such states are expected in these theories without any phenomenological input beyond 
PAMELA. Since Majorana fermions and real scalars cannot carry conserved charges, they 
cannot have diagonal couplings to gauge fields. As a consequence, there must be multiple 

1 See Ref. [1] for a broader discussion. 

2 This could be alleviated if the signal is produced from a nearby clump of dark matter [12]. 
3 The Sommerfeld effect [21] was first discussed in the context of dark matter by [22]. See also [23, 24] for 
other recent discussions. 
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Majorana fermions or real scalars in the theory, and the gauge couplings are off-diagonal. 
For the Sommerfeld effect to occur, these states must have splittings smaller than a typical 
WIMP kinetic energy, above which the enhancement is suppressed. Remarkably, with a GeV 
force carrier, the states are naturally generated from radiative effects in the ~ 1 MeV range 
[1]. Such splittings can lead to other signals with potentially observable ramifications. As 
we shall discuss, scatterings from the ground state into the excited state can provide simple 
explanations of the DAMA/LIBRA [30] and INTEGRAL [31] signals, using splittings of 
0(100 keV) and O(MeV), respectively. 

It has been previously noted [32] that excited states can have cosmologically interesting 
lifetimes and observable consequences, for instance on the CMB. These long-lived new states 
allow for a wide range of interesting phenomenology which we explore in this paper. In 
particular 

• Even with the large deexcitation cross sections that can arise with light mediators, it is 
possible to have significant relic populations of the excited states. 

• The lifetimes of these states are generically cosmologically long, although small cou- 
plings to the Z, or the presence of new, light states can shorten them dramatically. 

• There are strong direct detection constraints on the relic excited abundance from the 
unsuppressed down-scattering. If the abundance is near current limits, there may be a 
striking signal at upcoming Ge, Si and Ar experiments. 

• If the relic excited population is large, then the DAMA-associated splitting may arise 
from a metastable state transition into an excited state. That excited state then 
promptly decay into the true ground state via e + e~ emission. In this case, the "exciting 
dark matter" (XDM) proposal of [31] can accommodate lighter WIMPs. Such scenarios 
predict signals at direct detection experiments in the 800+ keVr range. 

Again, we emphasize that such excited states are automatically present in any theory with a 
new gauge interaction, and that the splittings of these excited states must be at most O(MeV), 
otherwise the Sommerfeld enhancement would be inoperative. As such, the motivation to 
consider inelastic up-scattering (as in [33]) or down-scattering from the excited state arises 
from the new theories related to PAMELA alone. Nonetheless, there are concrete reasons 
to consider specific scales for the excited states, associated with other possible dark matter 
anomalies, which we now review. 

1.1 DAM A and Inelastic Dark Matter 

The DAMA/LIBRA collaboration has recently reported a modulated signal at 8.2 a in their 
single-hit data, confirming their previous results from DAMA/Nal. The observed period and 
phase are that expected for WIMP dark matter. Unfortunately, no other experiment has con- 
firmed this signal, and if the signal is from elastic nuclear scattering, recent limits from CDMS 
[34, 35] and XENON10 [36] exclude this interpretation by over two orders of magnitude. 
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One proposal to reconcile these conflicting results is "inelastic dark matter" (iDM) [30]. 
In this model, dark matter transitions when scattering off of a nucleus into an excited state, 
with m x * — m x = 6 ~ lOOkeV. The kinematic difference causes heavy targets to be favored 
over light targets, a shift of the spectrum to higher energies (suppressing or eliminating low 
energy events), and enhanced modulation [30, 37]. Together, these weaken the bounds from 
other experiments sufficiently to allow the DAMA modulation to be consistent with the results 
of other experiments [38, 39, 40], even in light of the most recent results from XENON10, 
CDMS, ZEPLIN-III [41, 42], CRESST [43] and KIMS [44]. 

Models of iDM are simple to construct [30, 37]. Pseudo-Dirac fermions as well as complex 
scalars, with slight non-degeneracies between real and imaginary components naturally give 
rise to inelastic transitions. Majorana fermions charged under a non-Abelian gauge group 
[1] not only provide such phenomenology, but also explain the ~ MeV scale splitting if the 
mediator is ~ lOOMeV - lGeV [1, 45, 46]. 

However, it is essential for these models that the excited state is largely unpopulated. 
Deexcitations (down-scattering) are visible at all experiments, and for 0(1) populations of 
the excited state, would be significantly excluded. Smaller populations will provide for a 
remarkable signal at upcoming experiments. As we shall see in section 2, for the simplest 
models of iDM with a ~ GeV mediator, the excited state is cosmologically stable. This 
motivates us to consider metastable versions of iDM, where the excited state can decay to 
e + e~. 

1.2 INTEGRAL and Exciting Dark Matter 

The INTEGRAL collaboration [47] has reported a significant signal of 511 keV radiation from 
the galactic bulge region, approximately gaussian with a FWHM of 6°, with an additional 
subdominant component correlated with the disk [48]. Interpreting the signal as arising from 
the capture and annihilation of positrons and electrons in parapositronium, this requires 
3 x 10 42 e+/sec. 

The "exciting dark matter" (XDM) proposal of [31] postulated the existence of an excited 
state of dark matter x* with m x * — m x ~ MeV, which would decay via emission of e + e~. 
Such an idea takes advantage of the fact that the kinetic energy of a ~ 500 GeV WIMP in the 
galactic halo is approximately MeV. However, the large observed rate requires a large cross 
section as we now discuss. 

The total rate of excited WIMP states in the galactic center is simply 



where (a ex v) is the pairwise averaged excitation cross section, and rc is the maximum radius 
included in the INTEGRAL signal (approximately 1 kpc). We can calculate this by using 
the Einasto profile [49] with parameters set by the A-l run of the Aquarius simulation [50], 
which is one of the highest resolution simulations to date (although, importantly, we should 
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note it contains DM only). Doing so, one finds 
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We have not accounted for uncertainties in the local density here, which could add an addi- 
tional factor of a few to this estimate. This compares with the S-wave unitarity bound for a 
500 GeV WIMP in the halo, moving at v ~ 10~ 3 c, of approximately 3 x 10 _19 cm 3 s _1 . 

Unfortunately, not all of the particles in the halo can participate, as many pairs are below 
threshold. Taking a 1-D RMS velocity of v rms = 200 km/s one finds that approximately 1/20 
of the pairs are kinematically capable of scattering [31]. Thus, with these parameters, only 
cross sections which approximately saturate the unitarity bound will yield a sufficiently large 
rate. 

There are a few important comments to make here: first, the presence of a light mediator 
for this process can naturally allow a near saturation of the unitarity bound, even if the cou- 
plings, themselves are pertubrative [31]. This arises simply from multiple particle exchanges. 
Second, the above fraction assumes Maxwell-Boltzmann velocity distribution with a local ve- 
locity dispersion of v rms = 200 km/s, which is consistent with expectations arising from the 
galactic rotation. However, Maxwell-Boltzmann distributions underpredict the high-velocity 
component of the distributions when compared to simulations [50]. Moreover, v rms could 
change as one moves to the galactic center, and the results are exponentially sensitive to this 
value. Determining the properties of the WIMP velocity dispersion is very difficult in the 
galactic interior, because the potential is dominated by baryons. In fact, recent simulations 
including baryons show a rise of v rms as one moves to the center of the galaxy [51, 52, 53]. 
For v rms oc r -1 / 4 , as argued in [52], the velocity dispersion in the inner galaxy would be 
approximately twice as large, and the excitation process would be largely unsuppressed by 
the Boltzmann distribution. Hence, it is very plausible that even cross sections well below 
unitarity-satuation could explain the signal. Finally, recent measurements of the motion of 
the Milky Way [54] suggest that the rotation velocity is approximately 30 km/s higher than 
previously expected, additionally suggesting a mass 30% larger than previously thought, low- 
ering the needed cross section by a factor of 1.5 (from the density) and possibly significantly 
more (from the associated change in expectations for v rms , which naturally scales with v ro t). 
Thus, while we believe the original XDM proposal remains viable, it is interesting to consider 
the scenarios which allow the WIMP to be lighter 4 . 

2. Lifetimes and relic populations of excited states 

The thermal history of these dark matter models is somewhat different from a standard 
WIMP. Because of the new force carrier, which we denote with b, the WIMP stays in thermal 

4 Other authors [55], using an NFW profile, and using a smaller v rm s = 180km/s conclude differently. We 
note that the arguments that the velocity dispersions in the inner galaxy go down based upon [56] rely upon 
the assumption that the MW interior is dominated by dark matter, when it is, in fact, dominated by baryons. 
As we have noted, simulations involving baryons show the opposite trend [51, 52, 53]. 
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equilibrium via (a) XX bb and (b) be <-> 7e [31, 57, 32, 1, 45]. Chemical freezeout occurs 
when the former interaction (a) becomes inefficient, which happens at the usual T ~ m x /20. 
Kinetic decoupling from the SM bath occurs much later, however. The process xb — ► xb 
is efficient as long as b remains relativistic, and in chemical equilibrium with the standard 
model. Chemical equilibrium is maintained until T ~ mj,, because the reaction (b) is very 
efficient, with cross section a ~ a 2 e 2 m^ 2 [32]. As a consequence, kinetic decoupling can occur 
as late as 30 — 300MeV, or even later for lighter force mediators. After this point, the 
temperature of the dark sector scales as a~ 2 . Consequently, the DM can reach T ~ 100 keV 
before SM BBN occurs. 

In theories where the new force in the dark sector is identified with a spontaneously 
broken gauge symmetry, as in Refs. [1, 45], the situation has an added element. New degrees 
of freedom associated with the higgsing of the gauge symmetry are expected and can have 
long lifetimes themselves [1, 45, 46, 58, 59]. For example, a higgsino/gaugino-like state can 
have a long lifetime owing to its suppressed decay into a gravitino. Depending on the precise 
spectrum and nature of these states their relic abundance can span several orders of magnitude 
ilft ~ .1 — 10 5 . If such particles have lifetimes 10 4 sec, as they decay electromagnetically, 
their energy is harmlessly deposited into the photon bath [60]. However, in such a situation, 
the number density of particles off of which x* can down-scatter is many orders of magnitude 
higher. Therefore, the effects of these light scatterers on the population of the excited WIMP 
state can be substantial. 

2.1 The Relic Population of Excited WIMP states 

In the early universe, the populations of excited and ground states are equal, but as the 
universe cools, the number of excited states becomes Boltzmann suppressed, so long as de- 
excitation scatterings remain active. Indeed, for large (MeV) splittings, it is quite natural to 
have a small, but non-zero, relic fraction [32]. 

The degree to which the excited states remain present depends on a number of factors. 
First is simply the number of particles off of which an excited WIMP can collisionally deexcite. 
In some cases, x* can onr y deexcite via x*X* —> XX: m which case the number of scatterers 
falls off exponentially. In other cases, for instance with scalar or other interactions, x*X XX 
is also possible, and the rate for down-scattering is higher since the scatterer, x, maintains 
an unsuppressed number density. Finally, the presence of long-lived (r ~ 1 sec) light states, 
with possibly large population, offers additional scatterers against which x* can deexcite. 

A second important factor that determines the fractional abundance of excited states 
is the decoupling temperature of the WIMPs [32]. In particular, if the WIMPs kinetically 
decouple earlier, then the resulting relic fraction of excited states tends to be lower. This is 
so because after kinetic decoupling the temperature of the dark sector drops more rapidly 
than the photon temperature and therefore the deexcitation reaction remains in equilibrium 
longer. The upshot of all this, as we shall see, is that there is a wide range of possible excited 
fractions, from 0(1) to negligibly small. 
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To quantify this, we solve the Boltzmann equation in the limit where dn/dt = to obtain 
the temperature at freeze-out and hence the relic abundance. 

3H(z) = (a ex v)n s (z) (2.1) 

where H(z) is the Hubble parameter as a function of redshift, and ns(z) is the number 
density of downscattering targets. For self-scattering ns = n x * = n x e~ & / Tw where T\y is the 
temperature of the WIMPs (which can, in general, be different from the photon temperature). 
For more general interactions ns = n x , and for scatterings off of other long-lived particles X 
ns = nx- The fraction of excited states at freeze-out is then essentially just exp(— 8/Twf). 

The relic abundance of the excited population drops as the WIMP kinetic decoupling 
occurs earlier [32] (because the WIMP temperature scales as T = /Tkd, where Tkd is the 
temperature of kinetic decoupling). We thus consider two separate cases: first, we consider 
the case where T\y = T 7 , which should be thought of as a limiting case with the maximal 
relic excitation fraction. Because kinetic equilibrium is maintained generally by the thermal 
(chemical) equilibrium of the mediator b, this typically will not occur for splittings 5 ^ MeV. 
Alternatively, we consider Tkd = 1 GeV, which is a reasonable, if somewhat high value for 
Tkd- 

We further distinguish between the possibility that x*X* ~^ XX is the only deexcita- 
tion process, against the situation where x*X XX a l so contributes. In the former case, 
deexcitation terminates earlier as the excitated fraction drops exponentially. 

In Fig. 1 we plot the fractional abundance as a function of the thermal excitation cross- 
section (<J ex v) and the excitation energy, 5 for the case where x*X* ~^ XX ac ts alone and 
T\y = Tj. As can be seen from the plots, this fraction is 0(1), unless we consider very large 
thermal cross-sections, (<T ex v) > 10 - cm 3 s -1 or very large excitation energies 5 3> MeV. 
Notice that the fractional abundance is lower for a lighter WIMP. This is a consequence of 
the fact that a larger number density leads to a higher reaction rate and hence to a lower 
freeze-out temperature and further depletion. 

In contrast, in Fig. 2, we depict two alternative scenarios and the resulting fraction. On 
the left we consider a situation where x*X ~^ XX is present as well and therefore the density 
of scatterers does not deplete exponentially. On the right, we consider a situation where the 
decoupling temperature is T& ec = 1 GeV. Both effects tend to wash out the excited state 
population. If the two effects are combined, they may result in negligible abundance. 

If there are any light states in the spectrum which are long-lived and can participate in 
the deexcitation process, then depletion of the excited state can be even more complete. Since 
their number density can be considerably higher than that of DM, the light scatterers lead 
to very efficient depopulation of the excited state. As can be seen on the left pane of Fig. 2, 
when n s /n x 3> 1 the resulting fractional abundance can be extremely small. In section 3 we 
give a simple example which realizes such a scenario. 

In light of the above, it is clear that the fractional abundance of the excited state can 
range over many orders of magnitude and depends on the interactions and particle content 
present in the dark sector. It is therefore important to consider the different alternatives and 
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Figure 1: Contour plots for the fractional abundance of the excited state after freeze-out, plotted 
against n s /n x *{a ex v) (assuming the number of scatterers n s is exponentially suppressed by e _5 / T ), 
and the excitation gap, S. We consider a WIMP mass of 100 GeV (500 GeV) on the left (right). 



their respective signatures. If a large population of the excited state is present, then bounds 
from current direct detection experiment place restrictions on the size of the deexcitation 
cross-section and/or the deexcitation energy. It also calls for current experiments to search 
for such deexcitation events as we discuss in section 4. Finally, in the models where the 
iDM and XDM transitions occur from a metastable state with a sizeable relic fraction, as we 
discuss below, the XDM scenario to explain INTEGRAL is easier to realize since the energy 
threshold is considerably lower. 

2.2 The Lifetime of Excited States 

If the excited state can only relax down through a gauge-boson in its sector which is kinetically 
mixed with the SM, then the lifetime of such a state is longer than the age of the universe, 
unless the energy gap \S\ > 2m e i ectron . Above the electroweak scale the only mixing term 
allowed is, 

C D ~b^B^ (2.2) 

where B^ v {b^ v ) is the hypercharge (dark abelian) field strength. The mixing parameter, e, is 
constrained to be e < 10~ 3 [61]. This operator, below the electroweak scale, leads to mixing 
with the photon as expounded on in Refs.[57, 1, 45]. It also results in mixing of the dark 
sector with the Z° vector-boson which now couples to the dark abelian current. Below the 
Z° mass we should integrate it out to generate a coupling between the abelian gauge-boson 
and the neutrino current. That coupling will render the excited state unstable, and allow it 
to decay through an off-shell emission of two neutrinos. The lifetime for that process is given 
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Figure 2: On the left is a contour plot for the fractional abundance of the excited state after freeze-out, 
plotted against n s /n x (a ex v) (assuming the number of scatterers n s is not exponentially suppressed by 
e~ s / T ), and the excitation gap, 8. On the right we assume n s is exponentially suppressed by e~ s / T as 
in Fig. 1, but with kinetic decoupling at T^d = 1 GeV. In both cases the WIMP mass is fixed at 500 
GeV. 
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where m n , m p , Mz, and is the mass of the neutron, proton, Z° vector-boson, and dark 
vector-boson, respectively. T" 1 is the neutron lifetime, and we have taken the dark gauge 
coupling to be a d = 1/137. 

A faster decay mode is for the excited state to decay through an off-shell emission of 3 
photons. Similar to the Euler-Heisenberg Lagrangian, below the electron mass one generates 
a gauge-invariant interaction between the dark gauge-boson and 3 photons, 

C D (b^FapF^ + fb^F a/3 F*e\ (2.4) 

where the tilde denotes the dual field and 9w is the Weinberg angle. We estimate the lifetime 
by scaling the neutron lifetime and adding an extra 1/167T 2 to account for the additional 
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phase-space suppression, 

= —^(^-Y(—Y 

8 x 10 19 sec ^OOkeV J VlO 3 / 

where 5 is the effective available phase space (which is shared between three photons), and 9 
parametrizes any additional suppression present in the coupling beyond that introduced by 
the kinetic mixing coefficient, e. In the models presented below 9 can be naturally small in 
which case the resulting lifetime is extremely long. The combination 9tjm\ is bounded from 
above by direct detection experiments. The precise value depends strongly on the deexcitation 
energy 5, but for 5 ~ 100 keV it is 9ejm\ < 10~ 3 /Mf (see Fig. 8 below). 

If the excitation of the metastable state to the excited state in the GC is to explain the 
INTEGRAL signal via the XDM mechanism, the excited state must relax to the ground state 
before it traverses a distance larger than ~ kpc. This places a constraint on its lifetime, 

vt* ~ 1 kpc t* < 10 14 sec (2.6) 

where r* is the excited state's lifetime and we assumed v ~ 10 -3 . The excited state relaxes to 
an electron-positron pair via an off-shell dark gauge-boson emission, and its lifetime is almost 
identical to that of the neutron, 

r*~ n \al) \ m\ ) \m n -mj ~ 900 sec llO^J { <) 

We again use 9 to parametrize any additional suppression in the excited to ground state 
transition. Using Eq. (2.6) we can place a lower bound 9 > 10~ 5 from requiring the transition 
to happen within a kpc of the GC. We note that the #'s utilized in Eqs. (2.5) and (2.7) are 
logically disjoint as they pertain to different transitions. 

2.3 Other Decay Mechanisms for Excited States 

Before proceeding, it is important to consider the viable possibilities for a more rapid decay of 
the excited state. Essentially, the appearance of any light states with appreciable couplings to 
the WIMPs can lead to rapid decay. A few options are neutrinos (through direct Z-coupling), 
axions, or right-handed neutrinos coupled to light states in the dark sector. 

A very simple possibility is to consider the presence of some state with hypercharge 
(something akin to a Higgsino, for instance), which the dark matter particle mixes with by 
an amount £. The excited state in this case can decay via an offshell Z boson to vv. Even 
for £ ;$ 10~ 5 and 5 ~ 100 keV, the particle will decay quickly on cosmological timescales. 

Another possibility is decays to axions. Although they are not the principle dark matter 
candidate in these models, they still serve the important purpose of solving the strong CP 
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problem. If the splitting between excited and ground states occurs radiatively (or some 
component of it does), then there can be a decay into axions with lifetimes 



where S a is the contribution to the splitting which occurs spontaneously, and f a is the axion 
decay constant. Even for 5 a ~ 1 keV, the lifetime is cosmologically short. 

Finally, we can consider the possibility of new, light states in the dark sector that couple 
to right handed neutrinos. The dynamics in the dark sector leading to vector-boson masses of 
~ GeV often involves some light scalar fields which could couple to the DM states as well. A 
neutral state in the spectrum (such as discussed in [58, 59]) may also couple to right-handed 
neutrinos directly. Mixing between the different light states leads to an effective coupling 
between the DM states and right-handed neutrinos. If right-handed neutrinos are sufficiently 
light, then the DM excited state can relax into them on a cosmologically short timescale. 

3. Scenarios and Models for metastable iDM/XDM 

One alternative to having inelastic dark matter excite from the ground state into a 100 keV 
state is to imagine that there is a large population of WIMPs in a metastable state just below 
the e + e~ threshold of 1022 keV. In this case, the inelastic transition would be into a state 
that could efficiently decay into the true ground state 5 . Such a model would also allow the 
XDM scenario to be viable for lighter (~ 200 GeV) WIMPs. Although these particles would 
not explain the ATIC/PPB-BETS excesses, they would still be able to naturally explain 
PAMELA and the Haze. It is unclear what the limits on deexcitations from a ~ 800 — 900 
keV state would be, as they are outside the search windows of CDMS and XENON. An 
examination of these data could likely exclude a deexcitation rate with similar cross section 
to that for the iDM excitation. Still, because they arise from different gauge bosons, with 
a priori different mixings with the photon, the cross sections can differ by many orders of 
magnitude, as we shall show. 

To demonstrate this phenomenology, we consider here some simple field theories that 
result in a spectrum with a metastable state and an ample WIMP- WIMP inelastic cross 
section, but with suppressed deexcitation cross section. These models easily arise when the 
dark sector gauge bosons are not equally mixed with the SM photon. The first example 
we will provide illustrates the possibility of a metastable state with the necessary couplings 
to scatter inelastically off nuclei to an excited state which is ~ 100 keV above it as in the 
iDM scenario. The excited state then rapidly relaxes to the ground state through the same 
transition responsible for the INTEGRAL signal (Fig. 3). 

Although it is incompatible with the iDM scenario, we will also provide an intriguing 
second example. Here, we will exhibit a scenario with a metastable state which can scat- 
ter both up and down through the same gauge-boson. This allows for annihilations of the 

5 This same spectrum was recently employed by [62] . 
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metastable state into the excited state which is unsuppressed by kinematics. As a result, 
the annihilations in the center of the galaxy which lead to the INTEGRAL signal are not 
restricted to a small fraction of the WIMP distribution, but rather occur between any pair 
of metastable states. 

Finally, we briefly discuss models with light scatterers that can deplete the excited state's 
population as argued in section 2. 

3.1 Radiative Splittings and the iDM Scenario 

We begin by addressing a scenario in which both DAMA and INTEGRAL are explained by 
a ~ 100 keV inelastic scattering from a metastable state (5 « 900 keV) into a more rapidly 
decaying excited state (8 « 1 MeV). The downscattering process is generally suppressed, 
because the different gauge bosons naturally mix with different strengths to hypercharge. 

We begin with a model with SU(2) gauge symmetry under which DM transform as a 
triplet. The Lagrangian also contains two scalar adjoints of SU(2), <p an d <fi', 

£ D \M X X + V{<t>,<t>') (3.1) 

where w lxv (B^") is the dark SU(2) (SM hypercharge) field strength, and A ~ TeV is some 
high energy threshold. 

We assume that the scalar potential, V(4>, <f/) is such that <j) and <p' get a VEV and 
{<pz) 7^ and (0 2 ) 7^ 0. The gauge symmetry is broken completely and all 3 gauge-bosons 
become massive. At low energies, the fermions receive a radiative correction that splits their 
masses according to, 

5M 3 = ^(<P' 2 ) (3.2) 

where g d is the dark SU(2) gauge coupling, a d = gj/^, and we ignored corrections of order 
{(f)) /M. In the limit where {4>' 2 ) = (or (^3) = 0) we restore an unbroken U(l) with X3 (X2) 
a neutral state which does not receive a mass correction. With some tuning, it is possible to 
achieve a spectrum (Fig. 3, left panel) which contains all the appropriate scales for iDM and 
XDM. 

The higher dimensional operators in the Lagrangian of Eq. (3.1) will result in mixing of 
W3 and w 2 with the SM hypercharge. The mixing parameter is order ((/>)/ A ~ 10~ 3 — 10~ 4 
if the gauge bosons are at ~ GeV. As shown on the right of Fig. 3 this mixing generates 
the transitions necessary to explain both DAMA and INTEGRAL in the iDM and XDM 
scenarios, respectively. 
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Figure 3: On the left pane we show the splittings induced in the SU (2) triplet by radiative corrections, 
as well as all the relevant couplings. The mixing of the w\ gauge boson to the SM hypercharge is 
further suppressed compared with the mixing of W2 and W3 . This leads to the possibility of deexcitation 
processes that can be seen in direct detection experiments. The inelastic scattering off nuclei is depicted 
on the right. 

At higher order, w\ is also mixed with hypercharge at order {<j>) {<))') /A 2 ~ 10 — 10~ 8 . 
This coupling will allow for deexcitation transitions in direct detection experiments which we 
discuss in the next section. Note that in this case, 6 = (</>') /A < 10~ 3 and the lifetime in Eq. 
(2.5) is considerably longer than the age of the universe. 

3.2 Up-Down Transition for INTEGRAL 

In view of the possible significance of metastable states, another intriguing possibility arises, 
where metastable states scatter X1X1 ~^ X2Xo- This can yield a kinematically unsuppressed 
signal for INTEGRAL, although it does not provide us an obvious interpretation for DAM A. 
We consider again an SU{2) model, but take the DM fermions to be a massive Dirac triplet 
while the scalar sector consists of two real triplets, (j), 4>' as before, 



When (j) condenses, (^3) / 0, it breaks SU(2) down to U(l) and renders the appropriate 
gauge-bosons massive. It also generates a splitting in the fermionic sector between the neutral 
and charged components (charged with respect to the unbroken U(l)), as shown on the left 
in Fig. 4. Notice that the Yukawa coupling needs to be fairly small to allow for the MeV 
splittings needed for INTEGRAL. 

At this point, the model contains dangerous elastic transitions as well as a long range 
force (^3 is still massless). Also, the excited state x + cannot decay down to x~ an d generate 
the electron-positron pair needed for INTEGRAL. The condensation of (j) 1 resolves these 
problems. It splits the Dirac pairs x^ Majorana components with a mass difference of 
((/>2) 2 /A and suppresses the elastic transition. It also softly mixes x + with x~ by an amount 
(4>' 2 ) 2 1 '(2/(^3) A) to allow for a decay through UI3 (which both states couple to independently). 




(3.3) 
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Figure 4: The spectrum resulting from the Lagrangian in Eq. (3.3) is shown on the left pane 
together with the relevant couplings. The up-down scattering associated with INTEGRAL is shown 
on the right. Notice that dark charge is softly broken by the condensation of <j>' , hence allowing the 
seemingly charge violating transitions. 

3.3 Models with Additional Light Scatterers 

As we saw above, a wide range of models is possible. If the population of excited states is 
significant, we can have the INTEGRAL and iDM transitions occur from a metastable state 
(i.e., the inverted spectrum of subsection 3.1), or from pairwise up-down scattering. However, 
if we employ a spectrum such as considered in Ref. [33] or [1], where the lightest two states 
are split from one another by ~ 100 keV, the relic population of the excited states may be 
non-negligible, and it is important to consider the full range of how many of these states are 
present in the universe today. 

The reason for this is simple: as we shall describe, deexcitation scatterings (i.e., exother- 
mic transitions) are capable of occuring at any direct detection experiment. While transitions 
from ~ 300+ keV states are generally outside of the range studied by direct detection exper- 
iments, ~ 100 keV state down-scattering would be very visible, and is thus very constrained. 
As we shall see, with cross sections a n ~ 10~ 40 cm 2 , relic fractions larger than ~ 10~ 2 are 
strongly constrained by CDMS. 

Fractions of this size are roughly what is expected from deexcitations in the early universe 
through WIMP scattering. However, the presence of additional light states can diminish this 
fraction even further, by providing additional scatterers to keep the WIMP excited states 
in equilibrium until later times. Such states must satisfy three conditions: 1) they must be 
efficient scatterers; 2) they must live long enough to be present at the time of decoupling, but 
not too long as to spoil standard cosmological history; 3) they must have sufficiently large 
number density (compared with that of the DM states). 

The most obvious candidate to continue to depopulate WIMPs is the LSP of the new 
sector. For instance, the supersymmetric abelian models discussed in Refs. [58, 59] contain 
Higgsino-like states which mix with the gaugino associated with the abelian vector super- 
multiplet. In general, this results in mass eigenstates which are complete admixtures with 
elastic coupling to the dark force. Hence they form efficient scatterers if their lifetimes are 
long enough to allow them to be present at the time of decoupling. 



13 



One possibility is that this LSP is simply stable, for instance in a gravity-mediated 
scenario, or a different state in the sector which is stable for a symmetry reason. Such a 
possibility of course requires that this particle is not the dominant component of the dark 
matter, but that is easily the case, as we show below. Its number density can nevertheless be 
sufficiently high as it is much lighter than the DM states. 

On the other hand, LSP decay with a long lifetime is natural in gauge-mediated models. 
If the lightest fermion is heavier than any one of the dark gauge bosons then it can decay 
into it with the emission of a gravitino. The lifetime associated with this reaction is, 

/300 MeV\ 5 / VF V 

r ^~ lsec (,^H (joowj (3 - 4) 

where b should be understood as the lightest fermion, which has some component of the 
lightest dark gauge boson b, and F is the F-term vacuum expectation value that breaks 
supersymmetry. 

Finally, we must achieve a large number density of this light scatterer. If b is heavier 
than the lightest boson, then its relic abundance is determined by its annihilation into the 
lighter bosonic states. For example, in the zero-velocity limit, the annihilation into two dark 
gauge-bosons yields a cross-section [63] , 

di hh\ sin 4 ^ 3 ffd y i y 

where (3 = (1 — m^/m~), and sin <f> is the mixing between the LSP and the gaugino associated 
with the boson into which it is annihilating. In this case, f2 % 1 is natural (although could 
be much smaller if the only gauge boson in the dark sector is the one to which the DM 
dominantly couples). 

On the other hand, if b is lighter than any of the bosons, its annihilation can only proceed 
into SM particles. Hence, it is suppressed by the kinetic mixing parameter and is given by, 



ml \m b J m£ \m z z J 

where we have used e 2 /mf ~ l/m| which is normalized to the DAMA signal, neglected the 
small suppression coming from the muon mass insertion, and taken a d ~ a. Since the relic 
abundance of these particles is set by this cross section, we expect £1 ~ 10 2 — 10 5 . In this 
case, there is some tension in the lifetime of the LSP, with decays into a SM photon and a 
gravitino the only channel kinematically allowed, which is further suppressed by e 2 compared 
to Eq. (3.4). Either a slightly higher mass or lower SUSY breaking scale would be required 
to secure a decay which is cosmologically safe. 

While the vector-boson in this sector is very short lived and decays promptly into leptons, 
the scalar excitations may take longer to decay if they are lighter than the vector-boson (either 
through a four-body (h — > 2e + e _ ), or a loop-mediated two-body decay (h — ► e + e~)). They 
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can also satisfy the conditions of being efficient scatterers if additional Yukawa couplings with 
the DM states are present 6 . 

In either case, significant populations of additional scatterers are present, which then 
further depopulate the excitated WIMPs. As a consequence, a wide range of excited fractions 
are possible, and thus a wide range of rates possible at upcoming experiments. 

4. Signals of excited states at direct detection experiments 

The presence of a relic population of excited states can have dramatic consequences for direct 
detection experiments. Because the down-scattering process x*N — ► xiV is unsuppressed 
by kinematics, this process is visible at any target, including light targets such as Ge and 
Ar, where inelastic up-scattering is highly suppressed 7 . Still, because the kinematics are 
different, the spectrum of events is peaked at high energies, similar to the up-scattering case, 
and thus motivates a broadening of the search window for all experiments. 
We can categorize the interesting scenarios into a few specific cases: 

• The down-scattering can occur from a 5 ~ lOOkeV state which is invoked to explain the 
DAMA signal from inelastic up-scattering. In this case, CDMS places strong constraints 
on the relic excitation fraction of the excited state, but strong signals are possible at 
future upgrades, as well as at Ar experiments. 

• The down-scattering can occur from a 8 ~ 900 keV state which is invoked to explain 
the INTEGRAL signal. In this case, the signal is outside the CDMS range of 10 - 100 
keVr. Because the energy transfer is so large, the simple form-factor approximation 
is not a good one, and many of the scatterings would be into excited nuclear states, 
which would decay promptly with photons. Such events would be vetoed at existing 
experiments, but it is possible that a small number of coherent scatterings would still 
occur at high energies. 

• Because the PAMELA/WMAP/ATIC signals compel us to consider a new force, which 
is naturally accompanied by splittings in the ^ MeV range, we are motivated a broad 
range of scales, generically. It is possible that other excited states in this range exist 
with significant relic densities, and their down-scattering signals might be unconstrained 
by current WIMP search energy ranges. 

We calculate the scattering off of a nuclear target in the standard way [65] recognizing 
the effects of inelasticity [30]. For a target of mass m N and an energy recoil E R , the minimum 
WIMP velocity required to scatter is, 

Pmin = / n = 

6 The Yukawa coupling would have to be adequately small so not to introduce too large of splittings between 
the states. However, this small Yukawa could be compensated for by the very large number of Higgses. 

7 Similar down-scattering processes were considered in [64] with regards to light (sub-GeV) DM particles. 
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Figure 5: In the left pane we depict the recoil energy spectrum expected in CDMS for a deexcitation 
transition with \5\ — 100 keV. On the right we show the corresponding predicted number of counts 
observed in CDMS as a function of the deexcitation energy gap, \5\. The WIMP-nucleon cross-section 
x fractional abundance was taken to be Fa = 10~ 42 cm 2 . The WIMP mass is 100 GeV, 300 GeV, and 
500 GeV, from top to bottom. The horizontal line at 5.3 counts marks the 90% Poisson confidence 
upper limit on the expected number of signal events. 

where \i is the WIMP-nucleus reduced mass and <5 < is the energy gap associated with the 
deexcitation. When |<5| ^> 100 keV, scattering is dominated by events with very large energy 
recoil which are outside the experimental energy window. However, if \5\ ~ 100 keV, then one 
must worry about existing bounds from CDMS as shown in Fig. 5. Such bounds should be 
properly interpreted as a limit on the product of cross section and relic excitation fraction. 

Experiments with Argon as the nuclear target promise to reach much larger exposures 
than currently available by either CDMS or XENON10. Since Argon is a lighter element than 
Germanium, the usual iDM scenario is generally invisible in such experiments. However, as 
discussed above, a metastable state can leave its mark by scattering off the Argon nucleus 
and deexciting into the ground state. In Fig. 6 we plot the recoil energy spectrum in Argon 
for \5\ = 100 keV and several choices of the WIMP mass. Also shown in that figure is the 
variation in the spectrum as 5 is increased. The bulk of the signal is shifted to the right and 
might be missed altogether if the observational window is too restrictive. In order to contrast 
with other targets, in Fig. 7 we plot the spectrum resulting from a 900 keV deexcitation 
against argon, germanium, and xenon. The lighter target experiments perform better in the 
case of deexcitation. In Fig. 8 we present the discovery reach of argon and germanium in the 
\5\ — a plane. At very large |<5|, the experimental sensitivity is limited by the fact that most 
scattering events happen with a recoil energy which is outside the energy window (typically 
between — 150 keV). It is therefore important to broaden the search window to allow for 
nuclear recoils with higher energies. 

We note that at very large \5\, it is not clear that the form-factor really models the 
behavior of the scattering correctly, since the energy transfer is becoming comparable to 
nuclear excitation energies. However, the lowest lying excited state of °Ar is at 1460 keV 
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Figure 6: In the left pane we show the Argon spectrum for S = —100 keV with WIMP-nucleon cross- 
section a = 1Q- 42 cm 2 . The WIMP mass is 100 GeV, 300 GeV, and 500 GeV, from top to bottom. 
On the right, we fix the WIMP mass at 500 GeV and plot the resulting spectrum for \S\ — 100 GeV 
(solid), \S\ = 200 GeV (dashed), and \S\ — 300 GeV (dotted). As \S\ increases, the spectrum drifts 
outside the typical recoil energy observational window. 

(excited states of 36 Ar and 38 Ar are higher). Since DM states with splittings 5^1 MeV are 
cosmologically unstable, we do not expect excitations of this state and there should still be 
a clean signal of deexcitation in Argon experiments. This is also relevant for silicon targets, 
since the lowest lying excited state of 28 Si is at 1780 keV (that of 29 Si is only slightly lower). 
Therefore, clean deexcitations may also be observed with Si targets and in fact we can expect 
10's of such events with the current exposure of silicon at CDMS for a 500 GeV WIMP with 
a 10 -40 cm 2 scattering cross-section against nuclei. In contrast, future large exposures of Ge 
and Xe would be expected to yield significant rates of nuclear excitations, because of the 
many excited nuclear states which are accessible. Nevertheless, because these photons are 
monochromatic, they could be studied, rather than simply appearing as a background source. 

Another possible signature a population of excited states can result in is the indirect 
detection of energetic neutrinos coming from the capture and subsequent annihilation of 
WIMPs in the sun or the earth. In the inelastic case of iDM, capture is more difficult because 
the WIMPs cannot scatter as efficiently anymore. In particular, considering that iron is 
lighter than germanium, the bounds from CDMS imply that it is not possible to capture 
WIMPs in the earth in the iDM scenario. This is no longer true if an abundance of excited 
states is present since they can scatter against any target. The resulting neutrino flux may 
be detected in upcoming neutrino telescopes experiments. 

Finally, we should note that these exothermic reactions can result in electron recoils 8 . 
These would fail cuts in standard WIMP searches, but could possibly show up in careful 
analyses 9 . These electrons would have a spread of energies set by the range of binding 
energies of the atom in question, but would still be naturally narrow if the excited splitting 

8 See [67, 64] for a relevant discussion involving light DM particles. 
9 We thank Chris Stubbs for comments on this point. 
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Figure 7: The spectrum for 6 = -900 keV and a WIMP mass of 500 GeV with WIMP-nucleon 
cross-section er = 10 -40 cm 2 . The targets are argon (solid-blue), germanium (dotted-orange) , and 
xenon (dashed-red) . Notice that in the case of deexcitation, targets with lighter elements actually 
yield larger rates. 
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Figure 8: The cross-section per nucleon x the fraction of excited state in the halo, F, is plotted 
against the splitting, \8\. The WIMP mass is fixed at 200 (1000) GeV on the left (right) pane. The 
shaded area is excluded by CDMS and hence a fraction F ~ 10~ 2 — 10~ 3 is required if we assume a 
cross-section of a = 10 _4o cm 2 as required by DAMA. Discovery curves for the planned CDMS 4 th run 
and Argon target experiments are shown for exposures of 400 kg-day and 1000 kg-year respectively 
[66]. The dashed curves assume a recoil energy window between 10 — 100 keV for CDMS and between 
35 — 100 keV for Argon, while the solid curves show the large increase in sensitivity if a window of 
— 1000 keV is utilized. The upper curves correspond to CDMS, the lower to Argon. 

were sufficiently high. 
5. Conclusions 

A wide range of experimental results compel us to consider dark matter with ~ MeV excited 
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states. Explaining PAMELA with new gauge forces requires new states of dark matter, and 
the effectiveness of the Sommerfeld enhancement requires they have mass splittings not much 
greater than the kinetic energy of a typical WIMP in the halo. Beyond this, we have other 
significant motivations for excited states from completely different signals. In particular, the 
presence of a ~ 100 keV excited WIMP state can reconcile DAMA with other experiments 
through inelastic collisions with a nucleus, while a ~ 1 MeV excited state can produce the 
INTEGRAL signal from WIMP- WIMP collisions in the galactic center. While states with 
5 ^ 2m e are expected to be cosmologically unstable, those with smaller splittings can naturally 
have lifetimes longer than the age of the universe. However, it is important to keep in mind 
that simple alterations to the models can lead to more rapid decays. 

If the excited states are indeed long-lived then the interactions in the dark sector can 
yield a relic fraction of excited states spanning many orders of magnitude, from 0(1) to 
negligibly small. If such states exist, the downscattering signal for Ge, Xe, Si and Ar would 
be remarkable. S ~ 100 keV states would be visible at all experiments, while <5 ~ 1 MeV 
would be easily visible at argon and silicon experiments. For Ge and Xe, the scatterings would 
produce large numbers of nuclear excitations, with associated photons. Coherent nuclear 
scatterings would lie outside of presently considered energy ranges, motivating a study of 
higher (Er^ 100 keV) energies. 

Models with metastable excited states can change our interpretation of the INTEGRAL 
signal. For instance, if the DAMA transition is associated with an upscattering from a 
metastable state ~ 900 keV above the ground state to a state ~ 1 MeV above the ground state, 
then the INTEGRAL signal can be generated by WIMPs much lighter than 500 GeV. Such 
models naturally have an associated down-scattering signal as described above. Similarly, it is 
possible that the INTEGRAL signal is generated from collisions of metastable WIMPs where 
one is up-scattered and one is down-scattered. Such models also expand the parameter space 
for lighter mass XDM models, but do not obviously admit an interpretation of the DAMA 
data. 

Such scenarios show the continued richness of possible signals at upcoming direct de- 
tection experiments. It is important to control backgrounds over the largest energy range 
possible, and to consider signals which peak at high energies, in addition to conventional 
WIMP scatters which have an exponentially falling spectrum. The next few years will sig- 
nificantly probe the parameter space of excited WIMP states associated with these higher 
energy scattering events. 

Note added: While this work was ongoing, [62] appeared, which also considered the inverted 
spectrum for XDM discussed in section 3.1, although without connection to iDM or nuclear 
scattering. 
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